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Abstract 

We investigate the Lorentz structure of the confinement potential through 
a study of the meson spectrum using Salpeter's instantaneous approxima- 
tion to the Bethe-Salpeter equation. The equivalence between Salpeter's and 
a random-phase-approximation (RPA) equation enables one to employ the 
same techniques developed by Thouless, in his study of nuclear collective ex- 
citations, to test the stability of the solutions. The stablity analysis reveals 
the existence of imaginary eigenvalues for a confining potential that trans- 
forms as a Lorentz scalar. Moreover, we argue that the instability persists 
even for very large values of the constituent quark mass. In contrast, we find 
no evidence of imaginary eigenvalues for a timelike vector potential — even 
for very small values of the constituent mass. 

PACS number(s): 11. 10. St, 12.40. Qq 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) is believed to be the correct theory of the strong in- 
teractions. Still, the description of hadrons in terms of the underlying degrees of freedom 
of the theory, namely, quarks and gluons, remains a daunting task. Although fundamen- 
tal approaches such as lattice gauge theory continue to improve, it is unrealistic to forsee 
lattice calculations of the hadronic spectrum in the immediate future. Yet, it may still be 
fruitful to use lattice QCD to determine properties and parameters of a phenomenological 
model. Nonrelativistic quark models have been quite successful in this regard. For exam- 
ple, quarkonium, the bound-state spectrum of a quark-antiquark system, is well reproduced 
by a phenomenological potential consisting of the sum of a one-gluon exchange component 
and a confinement contribution. The quantitative description of meson masses, their static 
properties and decay rates count among the many successes of the model. 

Lattice simulations have provided conclusive evidence that QCD is a confining theory. 
For large quark-antiquark (qq) separations, lattice QCD predicts a linearly rising potential 
between heavy quarks However, heavy-quark systems place the strongest constraints 
on the short- distance part of the qq interaction and, thus, remain fairly insensitive to the 
details of the confinement potential. Thus, it is safe to conclude that the nonrelativistic 
quark model reveals little about the nature of quark confinement in the region where the 
model is expected to be valid. The light-quark spectrum is sensitive to the long-range 
part of the potential but the assumptions of the model are questionable in this case. Many 
attempts have been made to incorporate relativistic corrections into the nonrelativistic quark 
model For the most part these have consisted of a relativistic dispersion relation 

together with relativistic {p/M) corrections to the nonrelativistic potential. It is essential 
to stress, however, that these dynamical corrections can only be effected after the Lorentz 
structure of the potential has been established. At present, very little is known about the 
dynamical origin of the confinement potential and even less about its Lorentz transformation 
properties. Understanding the Lorentz character of the confinement potential — by means 
of a stability analysis of the instantaneous Bethe-Salpeter equation — is the main focus of 
the present work. 

The starting point for most relativistic studies of the meson spectrum is the covariant 
Bethe-Salpeter equation Unfortunately, there are difficulties associated with using the 
Bethe-Salpeter equation (BSE) for bound states. Chief among these is the appearance of 
a relative-time variable due to retardation effects. Three-dimensional reductions of the full 
BSE, which eliminate the relative-time variable, abound. While most of these reductions 
attempt to preserve fundamental physical principles, there is still no method of choice. Here, 
we work within the framework of Salpeter's instantaneous approximation to the Bethe- 
Salpeter equation 

We believe that the use of an instantaneous approximation is an appealing choice for 
the study of the Lorentz structure of the confinement potential. We recognize that, because 
very little is known about the dynamical origin of the confinement potential, introducing 
non- instantaneous corrections into the formalism is highly model dependent. We do not 
pretend to make any contribution to this subject. Instead, we simply adopt the position 
that if an ad-hoc modelling of retardation must be done, ignoring it altogether is probably 
the best choice. Still, we recognize that a choice for the Lorentz structure of the potential 
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must be made before it can be incorporated into a relativistic equation. We should also 
mention that because the long-range part of the potential is the main focus of the present 
work, the short-range, one-gluon exchange (OGE) component will be ignored. If, instead, 
the goal becomes the realistic description of the meson spectrum, one would be forced to 
include the OGE piece. Note that since the OGE contribution has a clear dynamical origin, 
incorporating retardation effects now becomes important and largely model independent. In 
this case, the use of an instantaneous approximation might be harder to justify. 

The Lorentz structure of the confinement potential is usually assumed to be scalar 0. 
Yet, several authors have suggested that the hypothesis of scalar confinement is not correct. 
Lagae ^ , and independently, Gara, Durand, and Durand examined several models within 
the framework of the instantaneous Bethe-Salpeter equation and concluded that the confin- 
ing potential was not a scalar. Archvadze, Chachkhunashvili, and Kopaleishvili examined the 
light-mass spectrum and concluded that no stable solutions exist for scalar confinement 0. 
More recently, Miinz, Resag, Metsch, and Petry found that reliable solutions of Salpeter's 
equation for scalar confinement exist only for large quark masses and weak confinement 
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In this paper we present a stability analysis of Salpeter's equation in the pseudoscalar 
channel assuming, both, scalar and timelike Lorentz structures for the confinement poten- 
tial. The stability analysis stems from the recognition that Salpeter's equation is identical 
in structure to a random-phase- approximation (RPA) equation familiar from the study of 
nuclear collective excitations. This stability analysis reveals the existence of imaginary 
eigenvalues for the case of scalar confinement. Moreover, an analytic study suggests that 
the instability should persist even in the case of a very large constituent mass. In contrast, 
no instability is observed for timelike confinement. 

We have organized the paper as follows. In Sec. |T| we review some general properties 
of the instantaneous approximation to the Bethe-Salpeter equation. Special emphasis is 
placed on the RPA structure of the equation. It is in this context that we make contact 
with previous stability analyses of the nuclear many-body problem. In Sec. |Tl| we discuss 
our method of solution and present results for the mass spectrum of pseudoscalar states for 
a variety of quark masses. Finally, Sec. |V| contains our conclusions and directions for future 
work. 



II. FORMALISM 

Here, we briefly outline a derivation of Salpeter's equation, using Greens's function 
methods, for a fermion-antifermion pair of equal mass (the extension to different masses 
is straightforward). In addition, we show the equivalence of Salpeter's equation to an RPA 
equation. 



A. The Salpeter equation 



The starting point for our derivation is the four-point Green's function, or two-body 
propagator, defined by 



(1) 



3 



where ipa{x) are fermion fields in the Heisenberg representation, a,(3,X,a are Dirac spinor 
indices, and \l/o represents the exact vacuum wave function. The four point function describes 
the propagation of two fermions (FF), or a fermion- ant if ermion (FF) pair, and contains 
all scattering and bound-state information. Restricting the relative time variables in the 
Green's function, enables one to select out, either, FF or FF modes from the propagator. 

The two-body propagator, and ultimately the bound-state spectrum, will be generated 
as a solution to the instantaneous Bethe-Salpeter equation in the ladder approximation. 
This is the central assumption of the model. Thus, in this approximation the irreducible 
kernel is given by 



V{Xi,X2) = 1/(xi,X2) . 



(2) 



In Salpeter's instantaneous approximation there are two remaining time variables (and, thus, 
one remaining relative time variable). By selecting them according to 

(3) 
(4) 



t' 



one isolates the propagation of the FF, as opposed to the FF, mode. In this limit, Salpeter's 
equation for the two-body propagator becomes, 

Ga/3;A<x(xi, xa; yi, ya; uj) = gJJS.;^^(xi, xs; yi, ya; uj) 

+ J rfzi rfz2GjJ|),.g^(xi,Z2;zi,y2;t^)V^^,,;^Y(zi,Z2)G5'/3;A,,(zi,X2;yi,Z2;c^ , (5) 

where G^^^ is the free two-body propagator and to is the energy variable conjugate to {t — t'). 
A study of the analytic structure of G, by means of a Lehmann representation, reveals that 
the propagator contains singularities at the position of the bound-state poles. Salpeter's 
equation can then be obtained by picking up the residue at the bound-state pole 

xfa(xi, y2) = J dzi dz2G[,°^l.g^(xi, Z2; zi, y2; E)V^n;e.n> {zi, Z2)x|^,,(zi, Z2) , (6) 
where the Salpeter amplitude has been defined by 

Xf,(xi,y2) = (^o|^a(xi)^.(y2)|^E) . (7) 
To interpret Xaa^ expand the fermion fields in terms of a free single-particle basis: 

M^) = E ([f^k.(x)],6.(k) + [vu^)]jm) , (8) 



where 6s (k) and ds(k) are second-quantized operators, and t/ks(x) and Vks(x) are free single- 
particle Dirac spinors (see Appendix 0). In this representation, Salpeter's amplitude be- 
comes 



Aacr ' 



xi,y2) 



E 

kiSi;k2S2 



f/k,s,(xi)] [f/k2S2(y2)] (vl/o|&.,(ki)6^k2)|v]/^) + 

^k,.,(xi)l \v^,s,{y2)] (^o|<(ki)42(k2)|^i?) + 

f/k,s,(xi)l IV^.sM] (vI/o|6.,(ki)42(k2)|^i.) + 

Iic,.,(xi)ljf/k2.2(y2)l {^o\dl{k^)bl{k,)\^E) 



(9) 



Since the free single-particle solutions are known, all dynamical information about the two- 
body system is contained in the four probability amplitudes (^{'^o\bsi0^i)bl20^2)\'^ e) , • • 
However, because of the instantaneous approximation assumed for the two-body kernel, two 
of these amplitudes vanish 



E 

kisi;k2S2 



leading, in turn, to the following form for Salpeter's relativistic amplitude 

t^kisi(xi)l [Vk2S2(y2)l 5,,,,(ki,k2) + 

L J a \- -i a 

^kisi(Xl) Uk2S2iy2) /^sis2(kl,k2) 

L J a L J o" 

Note that we have introduced the following definitions 

S,,,,(ki,k2) = (^'o|&si(kl)42(k2)|^'E) , 

D,,,,(ki,k2) = (vl/o|<(ki)6L(k2)|^^) . 



(10) 



(11) 



(12) 
(13) 



Finally, one obtains Salpeter's eigenvalue equation in the center of momentum frame (ki 
— k2 = k) by projecting out the two dynamical amplitudes from Eq. (0) 



r r/k' r 

[+E-2E,]B,,,,{k) = j [ (k;.i,S2|K++|k';4,4)5,,,,(k') 



+ (k; Si,S2|K+-|k';s;,4)D4,.(k') 



(14) 



r dV r 

[-E-2Ek]D,,sA^)= jTTT^Y. (k;si,.2|V^+1k';s;,4)5,;4(k') 

+ (k;si,S2|K++|k';s'i,4)D,;4(k') 

where the matrix elements of the potential have been defined as 

(k; si, S2|y++|k'; s[, 4) = [f/,,(k)] [t/.2(-k)] l^c(k - k') [t/,.^(k')] [Us',{-\^) 
(k; si, S2|\/+-|k'; s\, 4) = f/,,(k) t/,,(-k) l^c(k + k') K;(k')] [K^(-k')' 



(15) 

(16) 
(17) 



This eigenvalue equation is identical to the one derived in Ref. |TT] for the two-fermion case, 
with the two-body interaction (in the FF case) obtained from charge conjugation 



Vr. 



+V , for scalar, pseudoscalar, axial vector; 
—V , for vector , tensor. 



(18) 



The formal derivation of Salpeter's equation was carried out to the end of identifying 
its algebraic structure. Note that Salpeter's equation can be cast in the following compact 
matrix form: 
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(19) 



where the matrix elements of the "Hamiltonian" are given by 



(k;si,S2\H++\k';s[,s'^) = (k; Si, S2I V++|k'; s'^, 4) + 2Efc(27r)35(k - k')5,^,' 5,^/ , (20) 

J- 1 ^2 

(k;si,S2|i/+-|k';s;,4) = {k; Si, S2\V+'\k'; s[, s'^) . (21) 
The central result from the present section is the recognition that Salpeter's eigenvalue equa- 



tion, as given by Eq. (0), has the same algebraic structure as an RPA equation |[Tl|-[T3 
The equivalence between the two equations will enable us to use the stability analysis devel- 
oped by Thouless, in the context of nuclear collective excitations [0], to study the Lorentz 
structure of the confinement potential. 

It is also useful to study the dynamical content of the different contributions to Salpeter's 
equation. For example, the V~^^ term describes the direct scattering between the con- 
stituents. In the nonrelativistic limit, the iteration of this contribution to all orders results 
in Schrodinger's equation. In contrast, the V~^~ term arises from a genuinely relativistic 
effect, namely, the coupling between positive- and negative-energy states. In particular, this 
term is responsible for generating the double Z-graphs in time-ordered perturbation the- 
ory. Other contributions, such as single Z-graphs, do not appear in Salpeter's instantaneous 
approximation. 

We will be comparing our results in three different limits. As mentioned previously, 
Salpeter's approximation is obtained, from the full relativistic Bethe-Salpeter equation, in 
the limit of an instantaneous interaction. A simpler limit, given by the Breit equation, is ob- 
tained by neglecting the double Z-graphs from Salpeter's equation (this limit is known as the 
Tamm-Dancoff approximation in the context of nuclear collective excitations). Still, in this 
limit one retains the relativistic corrections to the potential and the relativistic dispersion 
relation. Moreover, one also incorporates the full Lorentz structure of the potential. Finally, 
the nonrelativistic Schrodinger limit is obtained by neglecting all relativistic corrections. 
Formally, it can be obtained by enforcing the following limits: 

i) Ek^ M + 



2M ' 

n) {k;s,,S2\V++\k';s[,s',)^V{k-k'), 
in) (k; Si, S2|l^^~|k'; s^, S2) = (k; Si, S2|^~^|k'; s^, Sj) = . 



B. Stability analysis 

Having identified the algebraic (RPA) structure of Salpeter's equation, we now employ 



the same formalism developed by Thouless in his study of nuclear collective excitations 
We perform the stability analysis by using confining potentials having both scalar and 
timelike-vector Lorentz structures 

V(r)=arl^'^^' for scalar; 

\ 7°72 J for timelike, ^ 
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where a is the string tension, and we note that both potentials become identical in the 
nonrelativistic limit. Salpeter's, and in general any RPA-like, equation can be rewritten 
as a Hermitian eigenvalue equation for the square of the energy [^,0. This implies that 
while the square of the energy is guaranteed to be real, the energy itself might not. The 
appearance of solutions having i?^ < signals, in the context of nuclear collective excitations, 
an instability of the ground state against the formation of particle-hole pairs — a collective 
mode with imaginary energy can build up indefinitely. Thouless has shown that the stability 
of the nuclear ground state depends on the Hermitian matrix 



(23) 




being positive-definite, i.e., all its eigenvalues must be positive |I^|-p^. This condition is 
equivalent to requiring that, both, the sum and difference matrices 



H+- 



(24) 
(25) 



be positive definite |T5[ . In this form, the stability condition of Salpeter's equation is reduced 
to finding the eigenvalues of the two Hermitian matrices and H~ . The existence of a 
single negative eigenvalue, of either or H~ , is, thus, sufficient to signal the instability. 
It is this criterion that we adopt here to test the stability of Salpeter's equation. 

We concentrate our analysis on the pseudoscalar (J'^ = 0~) channel. For this case, 
Salpeter's equation, for the reduced amplitudes b{k) = kB{k) and d{k) = kD{k), takes the 
following form 



rc 

i+E - 2Ek) b{k) = / 
Jo 



dk' 



I /o ^. 1 {k\V++\k')b{k') + {k\V+-\k')d{k': 
-E-2E,)d{k)= ^|(A;|V+-|A;XA;') + (fc|V^++|fcV(fc'; 



(26) 
(27) 



In spite of the simplicity of the angular momentum content of this channel, the matrix 
elements of the potential are complicated by relativistic corrections 



(k\V++\k') 
{k\V+-\k') 



Ek + M\ /Ek' + M 



2Ek J \ 2Ek' / 
Ek + M\ [Ey + 



2Ek 



2Eh.i 



l + ClCl]v,{k,k')T2Uk'V,{k,k') 

Cl + Cl]v,{k,k')±2CkCk'V^{k,k') 



(28) 
(29) 



where the upper(lower) sign in the above expressions should be used for scalar (timelike) 
confinement. We have also introduced the kinematical variable 



k 



Ek + M 



(30) 



to quantify the importance of relativity since, in the nonrelativistic limit {(f^, (f^i —>■ 0), these 
expressions reduce to 
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{k\V++\k') V^o(fc, k') , {k\V+-\k') . 



(31) 



Thus, the appearance of the £ = 1 component of the potential is a consequence of relativity. 
The angular-momentum components of the potential have been defined by 



Vc{k, k') = (47r)2 / dr jc{kr)V{r)jc{k'r) 



(32) 



with jc{x) = xjc{x) being the Ricatti-Bessel function. For confining potentials, the above 
integral is ill-defined. Hence, in examining confinement in momentum space we employ the 
following regularization for the spatial part of the potential |lT7|,p!8|: 



V{r) 



are ^ = a 



92 



-rjr 



dif 



The Fourier transform of the potential is now well behaved and is given by 



V^(k - k') 



drf 



Ana 



(k - k')2 + T]^ 



(33) 



(34) 



Evidently, we are interested in studying the stability of Salpeter's equation in the limit of 
rj 0. The stability analysis requires the explicit evaluation of and V~ . These are 
computed with the help of Eqs. (EHf) and p9|) 



V+{k,k') = {k\V++ + 
V-{k,k') = {k\V++ -v^ 



\k') = Voik,k') , 

\k') = v,{k,k')i{kX) 



(35) 
(36) 



where we have introduced relativistic "correction" factors, separately, for scalar and timelike 
confinement 



UKk') 
Uk,k') 



M2 



kk' Vi{k,k') 
EkEk'Vo{k,k') 
kk' Vi{k,k') 
EkEk' ' EkEk,Vo{k,k') 



EkEy 

M2 



+ 



(37) 
(38) 



Note that these correction factors are different because of an important relative minus sign. 
The expression for is remarkably simple — no vestige of relativistic effects remain. The 
only relativistic corrections to are, thus, kinematical and fully contained in the kinetic- 
energy operator Ek- The eigenvalue equation for is, then, a simple "nonrelativistic" 
Schrodinger equation with a relativistic kinetic-energy term: 



{Enr - 2Ek) hMR{k) 



dk' 
(2^ 



Voik,k')bNnik') . 



(39) 



Since the lowest eigenvalue of is at least as large as twice the constituent mass, we 
conclude that, if present, the instability must arise from the relativistic effects on H~. 

The relativistic corrections to V~ are contained in the expressions given in Eqs. ( P?]) 
and (|38|). In addition to simple kinematical factors, these terms contain important dynam- 
ical corrections which depend on the particular details of the potential. In general, these 
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expressions are complicated. Yet, there is one limit for which they are simple and illumi- 
nating. This is the k' = k limit for which the singularity structure of the potential — and, 
thus, the dynamical realization of confinement — become manifest. Separately, Vq and Vi 
are singular in the 7] ^ limit: 



Vo{k, k' = k) 
Vi{k, k' = k) 



SttV 
SttV 

^2 



[1 + 0{r^'/e)] 



V 1 V 



2k^ 



Ak 



(40) 
(41) 



Yet, the leading singularities cancel in forming the ratio 
Vi{k,k' = k) 



lim 



lim 



rf ( rf 
^ ~ 2P I 4P 



ri^oVQ{k,k' = k) v->o 
This result generates the following simple corrections to the nonrelativistic potential 



(42) 



Uk,k) 
Uk,k) 



M^-k^ 
M^ + k^ 

M^ + k"^ 



+1 
-1 



if k/M 
if k/M 



0; 
oo, 



(43) 
(44) 



Note that in spite of their simplicity, these expressions dictate how the singularity structure 
of the potential is modified by relativistic effects. 

For timelike confinement, the singularity structure of V~ remains unchanged from its 
nonrelativistic value. In particular, for k' = k we obtain [as in Eq. (001)] 



V-{k,k) = Vo{k,k)Uk,k) 



[1 + oiv'/e)] 



(45) 



Away from the {k' = k) singularity the potential receives corrections from relativity. While 
these corrections are quantitatively significant, they do not lead to any important changes 
in the qualitative behavior of the potential. These small changes have been documented 
in Table |, and also in Fig. |I], and support the assertion that a confining potential having 
a timelike Lorentz structure does not generate an instability. The picture is, however, 
drastically different for scalar confinement (see Fig. 0). In this case the singularity structure 
of the potential becomes 



V~{k,k) = Vo{k,k)Uk,k) 



■[l + Oiriyk' 



(46) 



where we have introduced a momentum- dependent string tension 



a{k' 



a 



—cr , 



if k/M 
if k/M 



0; 
oo. 



(47) 



The presence of this momentum-dependent string tension implies that, for momenta larger 
than the constituent quark mass, the potential shifts from a "rising" (a > 0) into a "sliding" 
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{a < 0) regime. Moreover, this behavior is unavoidable — provided one can mix in large- 
enough momentum components. However, if one limits the analysis to a small number 
of basis states, thus effectively introducing a momentum cutoff, the sliding regime will be 
missed — and so will the instability — whenever the constituent mass exceeds the value of 
the cutoff. This fact is clearly displayed in Table For a quark mass of M = 0.9 GeV, it 
is only after including 25 basis states that one can uncover the instability. 

We conclude this section by considering a mixture of scalar and timelike confinement, 
i.e., the Lorentz structure of the potential is assumed to be 

r = a;7?72°+(l-a:)lil2, (48) 
where x denotes the fraction of timelike structure. In this case 

V+ik,k') = Vo{k,k') , (49) 
" M2 kk' Vi{k,k'y 



v-{k,k') 



+ i2x - 1) 



EkEk' ' ' EkEk'Vo{k,k') 



(50) 



As before, is insensitive to relativistic correction and remains positive definite. In 
contrast, H~ is positive definite only for x > 1/2. Hence, any mix of scalar and timelike 
Lorentz structures has stable solutions only for x in the interval 0.5 < a; < 1. This result is 
in agreement with the numerical evidence presented in Ref. 



III. RESULTS 

For one-boson exchange potentials, like those encountered in the nucleon-nucleon prob- 
lem, one can evaluate Eqs. (^) and ( P7| ) quite efficiently in momentum space using a Gauss 
quadrature scheme ||TT|J18[| . However, it has been recently suggested that this representa- 
tion is not appropriate for the case of a confining potential due to its highly singular 
structure in momentum space. Evidently, the confining potential is most easily treated in 
configuration space. In contrast, the relativistic kinetic energy operator is difficult to handle 
in configuration space but easily treated in momentum space. We can accommodate both 
of these requirements by expanding Salpeter's equations in a suitable basis. Here, we use 
the radial eigenfunctions of the nonrelativistic harmonic oscillator, Rni, to expand the two 
amplitudes Blsj and Dlsj in terms of unknown coefficients 

BLSj{k) = Y.BnLSjRnL{k) (51) 

n 

DLSj{k) = J2DnLSjRnL{k) . (52) 
n 

This procedure results in a matrix equation for the unknown coefficients Bnisj and Dnisj 



which can be diagonalized using the method developed by UUah and Rowe |jT6[. Upon 
diagonalization one obtains the (previously) unknown coefficients from which one then can 
reconstruct the two amplitudes Blsj and Dlsj, and, ultimately, the Salpeter wave function 

We have made an effort to customize our numerical code so that it can be readily com- 
pared with previous results obtained from various approximations to the Bethe-Salpeter 
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equation |T9| , pO| . Specifically, mass spectra obtained from the Salpeter code have been 
checked against nonrelativistic (Schrodinger) results and in the so-called "spinless-Salpeter" 
limit where a relativistic dispersion relation has been employed. Comparisons have also 
been made to the results obtained by Long [jl2|, and Spence and Vary [0, in the Breit limit 
where the coupling to the negative-energy states is ignored {V^^ = 0). These comparisons 
were effected using a static potential consisting of a OGE piece plus a scalar confinement. 
We note that since these approximations neglect the coupling to negative energy states, they 
are not suitable to uncover the instability. Our results are also in agreement with a recent 
numerical study of the meson spectrum using the full Bethe-Salpeter equation (with a static 
interaction) in which it was found that the smallest possible eigenvalue decreases until it 
becomes imaginary as one increases the dimension of the basis Jl^ . 



We now proceed to show results for the meson spectra using a variety of constituent 
quark masses. Since the main goal of this work is the study of the Lorentz structure of the 
confinement potential, rather than a fit to the experimental spectrum, we ignore the one- 
gluon-exchange component of the potential. In Fig. ^we show the square of the ground-state 
energy in the pseudoscalar {J'^ = 0~) channel for a string-tension of a = 0.29 GeV^ and a 
constituent mass of M = 300 MeV. Results are presented as a function of the oscillator pa- 
rameter j3 for various sizes of the basis. We have indicated with "pluses" the results obtained 
by assuming a confinement potential having a timelike Lorentz structure. In this case, the 
development of a well-defined plateau signals the stability of the solution. In contrast, the 
presence of an instability in the scalar case (shown with "diamonds") is manifested by the 
appearance of imaginary eigenvalues {E'^ < 0) — even after including only five basis states. 

The presence of imaginary eigenvalues is only one of many possible manifestations of an 
instability. In particular, at the onset of the instability (i.e., for E = 0) the solution is also 
characterized by having B^sj = ^D^sj and, thus, zero norm. In Fig. |^ we have plotted 
the two (reduced) pseudoscalar amplitudes h and as a function of the relative separation 
between the quarks using 13 = 0.85 GeV. This represents a value of the oscillator parameter 
for which the energy is about to become imaginary. The approximate validity of the b = —d 
relation is evident. In addition, the behavior of both amplitudes, as well as the relative 
importance of the different basis states, suggests that the solution is, indeed, mixing in basis 
states having large momentum components. This result is in striking contrast with the one 
presented in Fig. ^ for timelike confinement. Not only do both amplitudes display a very 
smooth behavior with r, but, from the relative size of the amplitudes we can conclude that 
the relativistic corrections are small down to values of the quark mass of M = 300 MeV. 

In Figs. I^J^, and || we have carried out the same analysis as above but for a larger value 
of the constituent mass, i.e., M = 900 MeV. The main purpose of Fig. ^ is to indicate how 
easily one can be mislead into believing that the solution for scalar confinement is stable. 
Limiting the calculation to 5 or 15 basis states results in a smooth and stable behavior for 
the energy. Yet, as soon as basis states having large-enough momentum components are 
mixed in, the instability becomes manifest. The striking difference displayed in Fig. ^ for 
the behavior of the amplitudes as a function of the number of basis states, serves to further 
reinforce our conclusion. Finally, Fig. ^ shows the stability of the timelike solution against 
an increase in the number of basis states. 
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IV. CONCLUSIONS 



We have investigated the Lorentz structure of the confinement potential by means of a 
stabihty analysis of the meson spectra using Salpeter's instantaneous approximation to the 
Bethe-Salpeter equation. Because the algebraic structure of Salpeter's equation is identical 
in form to an RPA equation, we benefited from the techniques developed by Thouless in his 
study of nuclear collective excitations. The analysis revealed an instability, manifested by the 
appearance of imaginary eigenvalues, for a confinement potential assumed to transform as 
a Lorentz scalar. Moreover, an analytic study suggests that this instability will persist even 
for large values of the constituent mass — provided one can mix in large-enough momentum 
components. We have explicitly demonstrated that if, instead, the calculation is limited to 
a small number of basis states, one can conclude — erroneously — that the stability will 
get restored for large-enough values of the quark mass. In essence, we have shown that 
for momentum components larger than the constituent mass, relativistic effects modify the 
singularity structure of the potential and give rise, effectively, to a negative string tension. 

In contrast to the scalar case, we have found no instability for timelike confinement. Fur- 
thermore, a study of the lowest mass pseudoscalar state suggests small relativistic corrections 
down to values of the quark mass of M = 300 MeV. 

Many issues remain to be addressed. For example, how sensitive are our results to the 
choice of an instantaneous approximation. The answer to this question is difficult because 
very little, if at all, is known about the dynamical origin of the confinement potential. Thus, 
throughout this work we have adopted the same static form for the confining potential 
revealed to us from lattice gauge studies. Still, the role of non-instantaneous contributions, 
such as retardation, remains an important open problem. 

A natural extension of this work could consist of a detailed study of the meson spectra 
using a two-body interaction having a timelike confining piece plus a one-gluon-exchange 
contribution. This Lorentz choice for the confinement is dictated, not only by the previous 
stability analysis, but, also, by the correct reproduction of the linear Regge trajectories []7|,p!0 



Unfortunately, a timelike confinement potential leads to a relativistic spin-orbit contribution 
of the "wrong" sign 0]. Indeed, the sign of the spin-orbit splitting seems to constitute our 
best evidence in support of scalar confinement 0. One could attempt to remedy this 
deficiency by introducing a mixture of scalar plus timelike confinement |Q] or some other 
imaginative solution [0. Regardless of the final outcome, it should be evident that the 
long-range part of the qq interaction will remain at the center of many future investigations. 
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APPENDIX A: SINGLE-PARTICLE EXPANSION OF FERMION FIELDS 



The concept of a particle-antiparticlc removal amplitude is basis-dependent; one is asking 
what is the probability of removing a pair having, for example, momenta ki and \s.2, from 
the system. To address this, the fermion-field operators can be expanded in terms of an 
(initially) unspecified basis: 



(Al) 



The second-quantized operators bi,dl respectively annihilate a particle or create an an- 
tiparticle in the corresponding single-particle state; they satisfy the usual canonical anti- 
commutation relations. U, V are solutions to a single-particle Dirac Hamiltonian; they are 
orthonormal 



(A2) 



(A3) 



d2a;?7/(x);7j (x) = J d%;V/(x)Vj(x) = 6ij, 
J d^xUj(:K.)Vj(x.) = J dW/(x)C/^(x) = 0, 

and satisfy the completeness relation 

Y: [Ui{^)Uj{y) + V,{^)vHy)] = 5(x - y)l. 



We make use of the eigenstates of the free Dirac Hamiltonian, here defined as 
U^si^) = e*''-^t/(k, s) ; Vk.(x) = e-*'^-^V(k, s) , 



where 



U{k,s) 



Vik,s) 



Ek + M' 
'2Eh 



2E. . 



Xs -I 



\ Ek + M J 



1 / a-k 

2 



Ek + M \Xs, 



(A4) 



(A5) 



are positive and negative energy plane- wave spinors, Xs = {~)'^~^^X-s, with Xs the conven- 
tional two-component Pauli spinors 



' ^-h 







(A6) 
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APPENDIX B: PARTIAL- WAVE DECOMPOSITION 



With interactions that conserve total angular momentum and parity, it becomes conve- 
nient to rewrite Eqs.(0) and (^) in an angular momentum basis in which the Pauli spins 
si,S2 of the two particles are coupled to form a total S, then coupling the relative angular 
momentum of the two particles L to 5 to form the total angular momentum J of the bound 
state. One can then solve the equation separately for each different total angular momen- 
tum J and parity vr = (—1)^'''^. We will illustrate the partial-wave decomposition of the 
two-body amplitude with the direct term V~^~^ [Eq. (p!?])]. Including Dirac indices, this term 
is 



{k■,s^,S2\V^'-\k'■,s[,s',) = [f/.,(k)] Jf/,,(-k)]^[\/(| k-k' 



where the Fourier transform of the potential is given by 



f/,.(k')] ,[^4(-k') 



CMc 



/3' ' 

(Bl) 



[B2) 



In order to construct states of good total angular momentum, we write the free two-body 
state as a direct product of Pauli spinors, i.e.. 



[Usmjj,,{-k)i 



Ek + M' 
2Eh 



( 



a 



1 

-a ■ k 



\s1S2) 



li 



C^p{k) ^(aO; /30|A0) ^Ak){a^ap)^Js,S2) , 

A 



(B3) 



where Capik) is defined as 



Co.p{k) = V47r(-1) 



E. + M 



2E, 



Uk)Uk) ; Uk) 



Ek+M 



if a = ; 
if a = 1. 



(B4) 



Equation ( P3|) can now be combined with the partial-wave expansion of the Salpeter ampli- 
tudes [Eqs. ( pTOl) ] to give 

J2[U,,{k)UUs,{-k)]^B,,sM= J2 J^2iLSjik)BLSj{k\CSJM), (B5) 

siS2 LSCSJM 

where we have defined 

-a/3 



m,LSjik) = C^p{k)J2{aO;[30m{CSJ\\[Y,{a^a^),\j\LSJ). 

A 



(B6) 



Salpeter's equations in momentum space are then written in the angular momentum basis 
as 

[+E-2E,]hLSj{k)= j^,Y.\ (k-LSJ\V+Ak'-L'S'J)hL's'j{k') 



+ (fc; LSJ\V+-\k'; L'S'J)dL'S'Ak')^ , (B7) 
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rOO rlh' ^ . ( 

[-E - 2E,] dLsAk) = ^ (A:; LSJ\V-+\k'; L'S'J)bL's'Ak') 

Jo [i-Ky^ fj^, I 



+ (fc; LSJ\V—\k'; L'S'J)dL's'Ak')^ , (B8) 



with h{k) = kB{k) and d{k) = kD{k). The partial- wave expansions of the Salpeter amph- 
tudes B and D are given by 



BLsAk) = {LMl;SMs\JM) j dkYl,,^{k)BsMsiy) , (B9) 

Ml Ms 

DLsAk) = Yl {LMl;SMs\JM) J dkYl^jA^)DsMsiy) , (BIO) 



Ml Ms 



(Bll) 



where we have used the following definitions: 

BsMsi.^) = Y{siS2\SMs)B,,,,{k) , (B12) 

DsMsi^) = J2i-'iS2\SMs)i-)'-''-'W_,,_,,ik) . (B13) 

SlS2 

The amplitudes Blsj and D^sj represent the probability amplitude of removing and adding 
a particle-antiparticle pair to or from the system, with relative momentum k; within the 
framework of Salpeter's approach, they contain all dynamical information about the nature 
of the two-body bound state. The Salpeter wave function x in the angular-momentum basis 
is then given by 

Xal(fc)]^^, = i:{m;CSjik)BLsAk) + i-f^fs;CsAk)DLsAk)] . (B14) 

For local interactions, the matrix elements of the potential are given by (a sum over 
greek indices is implicitly assumed, and a = 1 — a) 

{k;LSJ\V++\k';L'S'J) = {k-LSJ\V—\k'-L'S'J) 

= Y.(-^T^'niLsAk){S\\[Vc{k,k')]^^.^^^ , (B15) 

cs 

{k;LSJ\V+-\k';L'S'J) = {k; LSJ\V-+\k'; L'S'J) 

= E i-ir^'^'^fs;LsAk) [Vcik, fc')]„,;„.,. I \S)T2s%s'Ak') ■ (B16) 
cs 

The quantum numbers L, S range only over the values allowed by J'", while C,S can take 
on all values allowed by the coupling to J. For E > 0, the amplitudes b, d satisfy the RPA 
normalization condition [|TB| 

roo rlh r -I 

E/ 7&sKj(k)-dlsAk)]=l. (B17) 



LS 



(27r)3 
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FIGURES 



FIG. 1. The relativistic correction factor ^y{k,k') in the hmit of 77 for a quark mass of 
M=1.0 GeV. 

FIG. 2. The relativistic correction factor ^s{k,k') in the Hmit of — > for a quark mass of 
M=1.0 GeV. 

FIG. 3. Square of the bound-state mass vs the oscillator parameter (/?) for a quark mass of 
M=0.3 GeV and a string tension of cr=0.29 GeV^. Results are presented for 5 (dotted line), 15 
(dashed line), and 25 (solid line) basis states using scalar (diamonds) and timelike (pluses) Lorentz 
structures. 

FIG. 4. Salpeter's b{k) (solid line) and d{k) (dashed line) amplitudes for scalar confinement 
using a quark mass of M=0.3 GeV and a string tension of £7=0.29 GeV^ Results are presented 
for 25 basis states using an oscillator parameter of /3=0.85 GeV. 

FIG. 5. Salpeter's b{k) (solid line) and d{k) (dashed line) amplitudes for timelike confinement 
using a quark mass of M=0.3 GeV and string tension of cr=0.29 GeV^. Results are presented for 
25 basis states using an oscillator parameter of (3=0.70 GeV. 

FIG. 6. Square of the bound-state mass vs the oscillator parameter (/3) for a quark mass of 
M=0.9 GcV and a string tension of cr=0.29 GcV^. Results are presented for 5 (dotted line), 15 
(dashed line), and 25 (solid line) basis states using scalar (diamonds) and timelike (pluses) Lorentz 
structures. 

FIG. 7. Salpeter's b{k) (solid line) and d{k) (dashed line) amplitudes for scalar confinement 
using a quark mass of Af =0.9 GcV and string tension of (T=0.29 GeV^. Results are presented for 
15 and 25 basis states using an oscillator parameter of /3=0.30 GeV. 

FIG. 8. Salpeter's b{k) (solid line) and d{k) (dashed line) amplitudes for timelike confinement 
using a quark mass of M=0.9 GeV and string tension of £7=0.29 GeV^. Results are presented for 
15 and 25 basis states using an oscillator parameter of (3=0.30 GeV. 
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TABLES 



TABLE I. Square of the pseudoscala 


r bound-state mass 




and the lowest 


eigenvalue of 


the matrix 


as a function of the number of basis states for 


various values of the 


mass. The 


calculation were 


done assuming timelike confinement and an 


oscillator parameter of (3- 


=0.25 GeV. 


Ma'S'j [A/TpV1 
ividoo iivj-C V 1 


^ states 


hi [vre V J 




\+ rr~'QA7'i 
A [<jevj 




/\ 1 vjrc V J 


ouu 





A zLOQ 












iU 


A QfiO 

4.0DZ 




i.yiz 




9 97R 




1 K 

io 


4.O0O 




1 01 1 

i.yii 




2.277 




on 


A OKT 




i.yii 




9 977 




OK 

zo 






i.yii 




9 977 


(JUU 


5 


5.330 




2.195 




9 zL98 




10 


5.262 




2.177 




9 zll 7 




15 


5.256 




2.175 




Z.41d 




20 


5.255 




2.175 




2.416 




25 


5.255 




2.175 




2.416 


900 


5 


8.413 




2.844 




2.958 




10 


8.265 




2.814 




2.937 




15 


8.251 




2.811 




2.935 




20 


8.249 




2.810 




2.935 




25 


8.248 




2.810 




2.935 



TABLE II. Square of the pseudoscalar bound-state mass E^, and the lowest eigenvalue of 
the matrix as a function of the number of basis states for various values of the mass. The 
calculation were done assuming scalar confinement and an oscillator parameter of /3=0.25 GeV. 


Mass [MeV] 


# states 


E^ [GeV2] 


A+[GeV] 


A- [GeV] 


300 


5 


-0.799 


1.925 


-0.188 




10 


-11.898 


1.912 


-1.739 




15 


-26.945 


1.911 


-3.096 




20 


-44.022 


1.911 


-4.292 




25 


-62.387 


1.911 


-5.371 


500 


5 


2.864 


2.195 


1.027 




10 


-0.916 


2.177 


-0.137 




15 


-11.240 


2.175 


-1.301 




20 


-24.254 


2.175 


-2.372 




25 


-39.003 


2.175 


-3.359 


900 


5 


6.821 


2.844 


2.362 




10 


6.752 


2.814 


2.231 




15 


6.752 


2.811 


1.513 




20 


6.751 


2.810 


0.733 




25 


-0.465 


2.810 


-0.041 
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